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Abstract 


A theoretical procedure to determine the heat transfer characteristics of 
the interline region of an evaporating meniscus using the macroscopic optical 
and thermophysical properties of the system is outlined. The analysis is based 
on the premise that the interline transport processes are controlled by the 
London-van der Waals forces between condensed phases (solid and liquid) . The 
procedure, is used to compare the relative size of the interline heat sink of 
various systems using a constant heat flux model. This solution demonstrates 
the importance of the interline heat flow number, A h^ V \ which is evaluated 
for various systems. The heat transfer characteristics of the decane-steel 
system are numerically compared with those of the carbon tetrachloride-quartz 
system. In order to evaluate the theoretical models and obtain direction for 
their refinement, the length of the interline region and the stability of the 
meniscus should be experimentally determined as a function of the heat flux 


and the interline heat flow number. 
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Nomenclature 


dispersion constant [J] ; 


Hamaker constant [see Eq. (III.2)][J]; 


mean separation between close packed planes [m] 

-2 

constant in equation (II. 5) [m ]; 


■&2TT, 


Planck’s constant (J*s]; 


-2 -1 

heat-transfer coefficient [W*m *K ] ; 


latent heat of vaporization [J*kg ] ; 


constant in Eq. [III. 18] ; 


curvature [m j ; 


r -1 - 1 , 

thermal conductivity [W*m *K ] ; 


adsorption constant; 


molecular weight [kg*mol ] ; 


-2 -1 

mass flux [kg*m *s ]; 


refractive index; 


pressure [N*m 1 ; 


heat transferred [W* 


heat flux [W*m 1 ; 


universal gas constant [J*mol - K ] ; 


radius of meniscus [m] ; 


solid thickness [m] ; 


temperature [K] ; 


-2 -1 

overall heat transfer coefficient, [W*m *K ]; 


velocity [m*s ] ; 


3 -1 

molar volume [m *mol ] ; 


Greek Symbols 


see Eq. (III. 2) ; 


adsorbed per unit area [mol'm ] ; 


dimensionless heat-transfer coefficient; 


interfacial free energy, [J*m ]; 


film thickness [m] ; 


dielectric constant, £ = £' + i£"; 


dimensionless film thickness; 


contact angle, [degrees]; 
absolute viscosity [kg*m ^*s ]; 


i — i 


f I I I I I ( 

1X1 

Greek Symbols (continued) 

. . r 2 -1, 

V , kinematic viscosity [m ■ s ] ; 

E, , dimensionless film length coordinate; 

E, imaginary part of complex frequency [rad's j 

IT., defined by Eq. (HI. 13]; 

p, density [kg*m ^]; 

O, evaporation coefficient [dimensionless]; 

0 ), defined by Eq. (III. 4) [rad*s ; 

0), angular frequency [rad*s “] ; 


Subscripts 

E, 

e, 

2 ., 

Zv , 
NE, 
s, 

O, 

v. 


evaporating ; 
evaporative; 
liquid phase; 
liquid-vapor interface ; 
non- evaporating ; 
solid; 

evaluated at interline ; 
vapor phase; 


Superscripts 

d, dispersion; 

averaged; 
id , ideal ; 

*, differentiation with respect to E> 

/s , complex index of refraction; 


f 
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I. Introduction 

In the analysis of heat pipes it has been common practice to use a liquid 
transport number, N =■ y h /V , to evaluate the heat sink capability of the 
evaporating meniscus. In addition, the understanding of the evaporating 
meniscus can be enhanced by defining a "meniscus shape factor", R (K - K ), 
which is the product of the radius of the capillary and the change in average 
curvature of the evaporating meniscus (1) . This shape factor accounts for the 
fact that liquid flows as a result of a change in the meniscus shape. Although 
these two groups aid our understanding of the evaporating meniscus , they do 
not describe the transport processes occurring in the interline region (junc- 
tion of solid, evaporating meniscus and vapor) . The analysis of the evapor- 
ating meniscus is extended herein to the interline region by evaluating the 
effect of the London-van der Waals forces between the solid and liquid on the 
pressure in the liquid and on the vapor pressure of the liquid. This analysis 
is based on the premise that the interline transport processes are ultimately 
controlled by the dispersion forces resulting from the electromagnetic fields 
of the solid and liquid. This effect becomes important when the thickness of 
the liquid is less than 10 7 m. Since this region is relatively short in 
extent, it accounts for approximately 10% of the evaporation process. 

However, since this is the location where the meniscus "hangs on to the solid", 
the stability of the meniscus is probably determined in this region. In addi- 

f 

tion, the contact angle fixes the extent of the meniscus. Therefore, the 
region is of critical importance. 

In order to quantify the effect of the liquid-solid system properties on 
the interline heat transfer coefficient a solution for a constant heat flux 
model of the interline region is obtained in Section II. This solution 
demonstrates the importance of the group of properties A h^/v, which we 
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call the interline heat flow number. This replaces the liquid transport 
number, Yh_ /V , as the important transport grouo in the interline region 

fg 

where the dispersion constant, A , is more important than the surface 
tension, y . 

Assuming that bulk values can be used for the heat of vaporization and 
the kinematic viscosity, the interline heat flow number can be evaluated 
if values for A are obtainable. This problem is addressed below in Section III. 
The literature on London-van der Waals forces has been searched and procedures 
to determine the dispersion constant A are presented. Numerical values of 
the theoretical interline heat flow number and the interline heat sink 
are then determined and discussed for wetting systems. 

In Section IV, the theoretical interline heat sink is determined for 
various liquid-solid systems. These results demonstrate the use of a pro- 
cedure whereby the basic macroscopic optical (electromagnetic) properties 
of the liquid-solid system can be used to predict and compare the heat 
transfer characteristics of the evaporating meniscus in the interline 


region . 
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II. Constant Heat Flux Model for Interline Region 

In a previous paper (2) . a procedure to determine the theoretical heat 
transfer coefficient for the interline region (junction of vapor, adsorbed 
evaporating thin film and adsorbed non-evaporating thin film) of a wetting 
film was developed. Since the analysis was based on the assumption of a 
constant liquid-vapor interfacial temperature, the effect of the thermal con- 
ducitivity of the solid on the process was not included. Herein, an approx- 
imate procedure to determine the heat sink capability of the interline region 
is developed in which the effect of the solid resistance to heat transfer is 
included. The formulation uses a simple one-dimensional constant heat flux 
model. Although the general approach used in Ref. (2) is continued, the 
solution herein is directed towards a limiting case in which the solid con- 
ductance is controlling. Since the value of the ideal liquid-vapor inter- 
facial heat transfer coefficient is very large, solid conductEince is found 
to dominate the overall heat transfer coefficient in the interline region 
after a small increase in the film thickness. The formulation allows the 
relative effect of the overall heat transfer coefficient and the "surface 
transfer coefficient" to be evaluated. The results are useful for predicting 
and comparing the heat sink capability of various liquid-solid systems for 
use in compact heat exchangers such as the grooved evaporator and the inverted 
meniscus heat pipe. An effective heat transfer coefficient is defined for this 
purpose. For a review of the prior literature the readers are directed to the 
extensive literature cited in Ref. (2) . 
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II. 1 Development of the Equations 

The system to be studied consists of a thin, evaporating, adsorbed film 
of pure liquid on a horizontal solid substrate of thickness, S , (see Fig. II. 1). 
The system is not drawn to scale in Fig. (II. 1) since S>>x . The shape of the 
thin film is steady since the evaporating liquid is replenished by fluid flow 
towards the interline from a source outside of the analyzed region. The 
temperature of the solid surface at y = -S is a constant, T^ , which is 
greater than the vapor temperature, T^ . If the film is sufficiently thin, 
it is kept from evaporating by London-van der Waals dispersion forces between 
the solid and liquid. If a section of the film is slightly tapered, portions 
of the film may be sufficiently thick to allow evaporation. The junction 
between the non- evaporating and the evaporating thin film is the interline. 

Fluid flows in the negative "X" direction in the evaporating thin film as a 
result of a pressure gradient in the film resulting from the film thickness 
gradient. 

Using the procedure developed in Ref. (2) the evaporative heat flux is 


q = 


h_ m 
fg e 




(II. 1) 


For rhe present case, q is assumed to be a constant, q , and equal to the 
average macroscopic heat flux through the system: 


U (T - T ) 
s v 



V 



(II. 2) 


where, the average overall heat transfer coefficient, u , is defined as 


U = 



(II. 3) 


in which, h^ is the average liquid-vapor interfacial heat transfer coefficient 
and 6 is the average thickness of the thin film. This approximate macroscopic 


T Ev\/ 



FIGURE (II. 1) Interline Junction of vapor, evaporating thin film 

and non-evaporat inq thin film on solid of thickness, S, (not drawn to scale) 
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formulation is primarily for the limiting case of S/k 2 »h^ ^ . Eq. (II. 2) 
is made dimensionless by using the previously defined variables (2) : 
ri = <$/<5 0 ; ? 2 = G X 2 ; n* = dn/d£ . 


( * v * 

n_ \ = w_ 

" / " h fg AG 



(II- 4) 



and the interline 


thickness , 

I 



<5 Q , is given by 

3 / TvTl; 

_ , / x, v 

V M h (T. — T ) * 

v fg lv v 0 


(II. 6) 


Since 


h_ A G 

-Jg 


V 



(T, 


£v 


" Vo = 


id 


(II. 7) 


this group gives the ideal heat flux based on the kinetic theory of vapor flow 
that would occur at a liquid-vapor interface due to the temperature differences , 
(T^ - T v ) q . Using Eqs. (II. 4 - II. 7) gives 

/ * \ * — 


n_ 

T1 



(II- 8) 


The dimensionless average heat flux, y , can also be written as 


Y = 

't - T 

IS Vi 

!_ 

id T id id 

h lv 5 h Jlv . hv 

( T !v * T v) 

0 

k *" k 0 + tr 

s a h £ v 


(II- 9) 


which explicitly includes the sum of the three relative resistances to heat 
transfer in the system: solid/ideal liquid-vapor interface; liquid film/ideal 
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liquid-vapor interface; liquid-vapor interface/ideal liquid-vapor interface. 

in Ref. (2) , T s + T^ , ^ ^ h Jv /k £ + 0 and T 3 h £v /ll £v * For this 

limit with a variable heat flux, the heat flux ratio depends on the extent of 

— - 1-2 

the zone of evaporation, £ (^l) / Y - 0.5 (r) + n ~) • Solving Eq. (II. 8 ) for 

constant 7 with the boundary conditions: 


(II. 10a) 


dn 1 
d? *0 


(II. 10b) 


we obtain 


in (n) - + n 0 £ 


(11. 11) 


For the present case, the dimensionless contact angle is taken to be zero, 


f) « 0 , and the extent of the region of evaporation, £ , is 


2 in 71 


( 11 . 12 ) 


„ „ 0.5 2 h, A Jin n 

2 &n r| _ fg 


(11.13) 


The size of the interline heat sink, Q , with a film thickness, U Y 1 , is 


Q = q X = q 


id 2 y in n 


(11.14) 


Using the previous definitions, Eq. (11.14) becomes 


-10 5 

Q = (2 (i t A (T - T ) V in U) 

* fg s v 


(II.15a) 
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the solution for the limiting case of constant interfacial temperature (2 ) 
is 

Q = h AG 0 ' 5 v" 1 [1-0.5 (if 1 + rf 2 n Cm) (II. 15b) 

fg 


Eg. (II. 15a) can be written various ways in order to emphasize different 
characteristics of the system. Rewriting Eg. (II. 15a) as 


Q = q x = 


h fg A ^ ^ 2 In n 


(II. 15c) 


demonstrates that the length x of the evaporating thin film is a measurable 
guantity which is inversely proportional to the interline heat sink, Q . 

We can also use these results to define an effective overall heat transfer, 



(II. 15d) 


in which U - k/S . 
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XI. 2 Discussion of Theoretical Equations 

Equation (11.13) gives the length of the evaporating film for given 

2 h R 

values of the dimensionless thickness, ri , and the parameter — 

_ ^ q _ 

in which q = U (T - T ) . In the limit of negligible heat flux, (q = 0) , 

an infinite length of an adsorbed non-evaporating film of thickness 6^ is 

present. In Eq. (II. 3) , U = 0 because h^ =0. To compa.be the heat sink 

capability, Q, of various systems using Eq. (II. 15a) with equal dimensionless 

thickness, l")r the parameter, [2 U h A (T g - T^) /V] , is important. This . 

includes the usual dimensions of the system and the properties of the liquid 

and solid. In addition, it includes the dispersion constant. A, which 

accounts for the London-van der Waals dispersion force between the solid and 

liquid. Therefore, the relative effect of the classical overall heat transfer 

coefficient, 0, and the "surface transfer coeffcient", A/V , can be evaluated. 

This also defines an effective heat transfer coefficient, [2 U A h _ /V (T - T )]°" 

fg s v 

based on the film thickness range Tl > 1. Further evaluation of Eq. (11.15) 
suggests that the important group is A h /v which we name the "interline 
heat flow number" . This gives the variation of the heat sink capability of 
a thin film of thickness range T) > 1 as a function of the solid-liquid system 
for a fixed product of U At = k At/S. This important group is evaluated 
in Section III. 

It is not possible to increase the heat flux without bound as, suggested 
by the equations because the system would become unstable. This can be seen 
by rewriting Eq. (II. 9) as 

u (T s - V =1 «ilv •' Vo 4v (II - 16> 

For the case of constant liquid-vapor interfacial temperature (negligible 
solid resistance) the heat flux ratio, y , has a limiting value of 
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Y = 1 - 0.5 (tT 1 + rf 2 ) (11.17) 

For the constant heat flux case, y , would be less than this value. The 

limit for the value of (T^ v - T^) Q depends on the ability of the evaporating 
film profile to change and provide the necessary pressure gradient for flow 
which is finite. This stability limit is discussed in Ref. (2) . 

In Figure (II. 2), the dimensionless heat sink capability, y £ , is 

presented as a function of the dimensionless thickness, g , for the following 

four cases: y = 0.5; y = 0.25; y = 0.05; y = 1 - 0.5 (r| ^ + 0 2 ) . The 

value of the heat sink for the constant liquid-vapor interfacial temperature 
case starts below the constant heat flux case because of effect of the disper- 
sion force on vapor pressure. However, it surpases the constant heat flux 
case at a very small dimensionless thickness. It is interesting to note that 
the value of y £ at r) = 1.1 is approximately 0.25 y £ (r| = 10) for y =■ 0.25. 
Therefore, the value of £ at ri = 1.1 is approximately 0.25 £ (p - 10). This 
emphasizes the result that the thickness of the thin film increases slowly at 
first and then progressively faster. The thin film profiles for three values 
of y are also presented in Fig. (II. 2). Comparing these profiles demonstrates 
the result that the thin film slope predicted by the constant interfacial 
f.e:aperature model is considerably smaller than that predicted by the constant 

r 

heat flux model. Combining the results of the two limiting cases for a closer 
approximation of the real physical situation, we can say, qualitatively, that 
the initial portion of the evaporating thin film fallows the constant tempera- 
ture model until thermal resistance limits the heat flux to a constant value. 
This limit occurs at a very small thickness increase. Numerical examples of 
the interline heat sink are given later in Section IV. The interline heat 


flow number is evaluated in the next section. 
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III. Evaluation of the Interline Heat Flow Number : A h, V 1 

fg 

The objective of this project is to quantify the effect of the liquid- 
solid system properties on the interline heat transfer coefficient. In order 
to do this, an effective heat transfer coefficient has to be defined that 
includes both the solid conductance and a surface transfer coefficient for 
liquid flow, A/v. As described above, this has been accomplished by using a 
constant heat flux model. The resulting effective heat transfer coefficient is 


found to be [2 0 A h, /V iT] 

"<3 


0.5 


In this coefficient the Interline Heat Flow 


Number, A h^ /V, replaces the Liquid Transport Number, y h^/V, as the important 
transport group. That is, the dispersion constant. A, becomes important instead 
of the surface tension when the meniscus becomes sufficiently thin so that the 
solid effects the liquid properties. Having defined the heat transfer coeffi- 
cient in this manner, the problem is reduced to evaluating the dispersion con- 
stant A which accounts for the London-van der Waals forces between condensed 
phases . 
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III.l The Dispersion Constant, A 

For the purpose of this study we assume that the vapor pressure over an 
adsorbed film of liquid is a function of the temperature of the liquid and the 
net London-van der Waals force of attraction between the solid and liquid ( 2 ) . 
The force of attraction between the solid and liquid can result in a lowering 
of the liquid vapor pressure, P , relative to the normal thermodynamic vapor 
pressure, at the liquid temperature, T^ . This vapor pressure lowering 

can be related to the dispersion force by the following equation: 


R T £v 0 P v£v 
ion 


A 


V 


a 


' v&vt 


(III.l) 


where 5 is the thickness of the absorbed layer and A is the dispersion constant. 

— 3 

In tl. is case, - A /6 represents the relative force of attraction per unit area 
between the vapor and solid across a film of liquid, (A > 0) . The van der Waals 
force of attraction between any two atoms or molecules may be divided into 
dispersion, induction and orientation forces. Orientation forces arise only 
in the case of polar molecules. A polar molecule may also induce polarity in 
a nearby neutral molecule and this results in an attractive induction force 
between them. With non-polar molecules or atoms the non-polarity is time 
average. At very short intervals the molecules and atoms are found to possess 
finite fluctuating dipole moments. These instantaneous moments result in an 
attractive force between the molecules or atoms by inducing polarity in neighbor- 
ing atoms or molecules. The attractive force between condensed. .bodies (solid or 
liquid) resulting from these fluctuating electromagnetic fields can be described 
using the macroscopic optical properties of the materials and are called London- 
van der Waals dispersion forces. Except for the case of strongly polar molecules, 
the dispersion force generally dominates. Herein we will be mainly concerned 
with non-polar molecules. A good description of the development of the London- 
van der Waals force concept is given by Chu ( 3 ) . 
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111,2 Calcut&tioft of the Dispersion Constant, A 

Basically, thdl'Q are two methods of calculating the London-van der Waals 
dispersion force between condensed bodies , The microscopic approach starts 
with, the interaction, bstv-'een individual atoms or molecules and postulates their 
aaditivity so that the attractive force between the macroscopic bodies can be 
calculated by integration Qt'gr all atoms or molecules ( 4 ) . The result is a 
product of a purely geometrical term and the Hamaker constant, A^^, which 
depends on the materials involved. The second method, the macroscopic approach 
develop^ by 2ii£ shitz and co-workers ( 5 ) , uses the macroscopic optical proper- 
tied Qf the interacting macroscopic bodies and calculates the London-van der Waals 
dispersion force from the imaginary part of the complex dielectric constants. 

It appears to be generally accepted in the literature that the macroscopic 
approach is physically more satisfying than the microscopic approach. Presently 
we are concerned with a first order approximation of the adsorption effect on 
heat transfer. Therefore it is advantageous to use the results of both calcula- 
tional methods in order to obtain maximum use of the previous literature . 

Where the results are different, preference will be given to the macroscopic 
approach . 

When the molecules are an appreciable distance apart the fluctuating dipoles 
are no longer in phase since the time required for the electromagnetic field of 
one instantaneous dipole to reach a neighboring molecule and return "may be com- 
parable to the fluctuating field itself. In this case, the dependency of the 

-3 -4 

force on distance changes from 5 to 6 and the interaction is known as the 
retarded London-van der Waals interaction. The transition is gradual. The 
case of small 6 is of greater interest herein and we restrict ourselves to 


the non- retarded attraction. 

A good physical description of the microscopic approach is given by 



15 


Sheludko ( 6 ) . Using the microscopic procedure the equation for the temperature 
independent dispersion constant is 


— _ ^&s _ . fi s&v 

6 it ~ Sir 


(III. 2) 


T 3 

W 
— 2 
v„ 


* 2 h t 

Als v £ V S 


where v, 1 is the number of atoms or molecules per unit volume and 3. . is the 
i ID 

constant in London ' s equation for the attraction between two atoms . If the 
force of attraction between the liquid and solid (represented by ) is 
greater than the force of attraction between two liquid molecules (represented 
by the liquid wets the solid and a decrease in the vapor pressure occurs. 

Using the macroscopic procedure, an equivalent Hamaker constant, A^^, 
can be obtained from the Lifshitz-van der Waals constant, ft : 


3 ft to 


(III. 3) 


Using this definition to relate both approaches through an equivalent Hamaker 
constant, the two procedures can be unified and compared. 

For the purpose of this report the Lifshitz-van der Waals constant can be 
calculated using the following approximate formula ( 5 ) : 


[^(U) - e s (i5n t£j,(iS) - 1] 

+ e s (i£H Ce^dC) + U 


(hi. 4) 


where £,£,£„ are the frequency dependent dielectric constants of the solid, 
s v a 

vapor (£ =1), and liquid evaluated along the imaginary frequency axis, i£ . 

(The assumptions associated with deriving this equation from a more general 
formulation are presented in Reference 5.) In turn these imaginary frequency 
dependent dielectric constants can be calculated from the imaginary part of the 
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complex dielectric constant which can be experimentally measured. The dielectric 
constant, £ (W) , is a complex function of the frequency, 0). 

£(W) = £'(W) + i £" (U» (III. 5) 

The function £ (to) is related to the complex index of refraction , n , the refrac- 
tive index, n , and the absorption constant, k, by 

/e (u) = n = n(tO) + i k((0) (III. 6) 

£' (to) = n 2 (oo) - k 2 ( 0 J) (III. 7) 

e"(u)) = 2n(to) k (to) (ill. 8) 


For the purpose of analysis, it is useful to consider the frequency, to , to be 
a complex variable: 

(0 = 0)' + i£ (III. 9) 


As presented above , the dispersion force depends only on the value of the 
dielectric constant on the imaginary frequency axis, E (iq) . The function is 
always positive and is related to the dissipation of energy in an electromag- 
netic wave propagated in the medium. £(i£) is a real quantity which decreases 

2 

monotonically from £ q (the electrostatic dielectric constant) £ q = n >1 for 
dielectrics or from + 00 for ideal metals at u = iO to 1 at 10 = i“ for all mater- 
ials ( 7) . without the need for a specific model, the Kramers-Kronig relations 
can be used to calculate E(i£) from the imaginary part of the dielectric con- 
stant, e"(w), ( 7). 


e(i£) 


1 




X £" (X) 

•2 > 2 . 
(x: +■£ ) 


dX 


(III. 10) 


where £ and X are angular frequencies of the electromagnetic field. Therefore, 
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e(i£) can be obtained directly from experimental data on the refractive index 
and adsorption constant. However, due to the lack of sufficient data for the 
frequency dependent optical constants, semi-empirical models are usually used. 
On the other hand, due to the basic importance of optical data in many fields, 
there is a wealth of procedures and models available in the literature. In 
addition, the volume of available data is also rapidly expanding. . 
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I II. 3 Contact Angle 

For the purpose of completeness , it is desirable to address the related 
topics of contact angle and surface tension. The Young-Dupre equation for the 
contact angle between a liquid, & , and a solid surface, s, is (8 ) 


Y H 005 6 " Y sv - Y £s 


(III. 11) 


cl cl 

In general, v = y + Z y , in which y is the contribution to surface tension 
11 m 1 m 1 

resulting from the dispersion forces and E y are the contributions due to 

m m 

other interactions (dipole-dipole, metallic, etc.). The surface free energy of 
the solid-vapor interface, Y sv t is related to the surface free energy of the 
uncontaminated surface by [papers in (9) discuss this material] 


y = y - tt 

1 sv 1 s e 


(III. 12) 


The difference between these two surface free energies is the decrease in the 

surface free energy of the solid surface, tt , resulting from adsorption 

0 


7T = RT 
e 


r 

1 


V 


r d in P 


v 


(III. 13) 


If the interaction between the solid and liquid is by dispersion forces only, 
the solid-liquid interfacial free energy can be represented by 


d d 


= + Y s " 2 \PS, Y s 


(III. 14) 


Combining these equations with tt ~ 0 gives 

0 


Y J. (1 + cos 0) = 2 Yg 


(III. 15) 


Finally, the expression for the contact angle becomes 


cos 0=2 



(III. 16) 
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In general, the contact angle can be related to the Hamaker constant by using 


Y £ K U A U 


(III. 17) 


in which 


is a constant related to the molecular dimensions. 


Using Eqs. 


(III. 17) and (III. 16) for a system in which y^ = y^ we obtaih 

/“ A K 

cos 0 = 2;/ ~ -1 (III. 18) 

V 

If A K > A K , 0 = 0, and the liquid forms a stable wetting film on 
ss ss ££ ££ 

the solid surface. These assumptions apply to the simple system of a non- 
polar fluid interacting with a surface through dispersion forces.- e.g. , an 
n-alkane wetting a smooth vitreous silica surface. This was confirmed in 
the experimental studies of Blake (10) and Ingram (11) on the wetting of 
a-alumina and silica by n-alkanes . 

A much more complicated fluid is water. Even though gold has a very high 
energy surface, it has been demonstrated that water does not wet un contaminated 
gold (12). After considerable research, the following explanation has been 

deveoped in the literature (12, 13). Since water is polar, y^ ^ y^ - Instead, 

2 d ^ d 2 

y - 0.0728 J/m , and y = 0.0218 J/m“. In addition, for gold, y = 0.121 J/m . 

Using Eq. (III. 16), the theoretical contact angle is 65.5° which has been confirmed 
experimentally (12). Therefore, the polar nature of water and .its attraction 
for itself prevents water from wetting uncontaminated gold. A further complica- 
tion occurs when the gold surface is oxidized. In this cause, the surface 
becomes hydrophilic because there is a larger dipole-dipole force of interaction 
between the metal oxide and the water. The above discussion concerning water 
has been included to emphasize the increased complexity introduced by the 
presence of strong permanent dipoles. A further complication can occur in the 
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case of "autophobic liquids", which are liquids that cannot spread on adsorbed 
oriented monolayers of their own vapor (14) . In this case, the adsorbed oriented 
molecules form a film whose critical surface tension to wetting is less than the 
surface tension of the liquid itself. In addition to these effects it is always 
possible to have a contaminated surface. For example, a thin layer of adsorbed 
water vapor can drastically alter the wetting properties of a high surface energy 
material for hydrocarbons. 

In summary we can say that, although the interline region is very compli- 
cated, sufficient theoretical models are available to describe the interline 
region. However, it is difficult to know the exact experimental conditions 
at the interline. Therefore, unless an effort is made to fully characterize 
an experimental heat transfer surface inconsistencies in the heat transfer 
and contact angle data will be prevalent. 
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III. 4 Representative Values of [£(i£) -11 

Equation (III. 4) demonstrates that a thin Wetting film of adsorbed liquid 

on a solid surface will be stable if the dielectric constant of the solid along 

the imaginary frequency axis is greater than that of the liquid, E s U.£) > e^(i^). 

Therefore it is instructive to compare the values of E (i£) for various materials. 

In Fig. (III.l), the &og [e(i£) -'ll is presented as a function of the energy 

associated with an electromagnetic wave in electron volts, £og 10 (h£) • In gen- 

2 

eral, E(i£) decreases monotonically from £ q = ft >1 (for dilectries) or 
+ co (for ideal metals) at w = iO to 1 at u = i® . The source of the individual 
values are given in Appendix A . 

For non-polar dielectric materials the static dielectric constant, E q , 

2 

can be obtained using the Maxwell relation , £ = n , in which the index of 

o 

refraction, n , can be obtained in the visible range of frequencies. It is a 
measure of the induced polarization. The n-alkanes and diatomic molecules are 
representative of these materials. Curves for pentane and helium are presented 
in Fig. (III.l). For comparison, limiting values of nitrogen, oxygen, methane 
and carbon-tetraf luoride are also given. 

Polar molecules have in addition to the induced polarization an orientation 
polarization. As a result the dielectric constant is not constant in the fre- 
cmency range below the visible. The alcohols and water are representative of 
this group. The curve for water in Fig. (III.l) stops at h£ = 1. At £ = 0, 
the dielectric constant for water has a value of £ = 81. At h^ = 3.16, 

£ Z 1.75. At most frequencies the dielectric constant for pentane is less 
than that for water. Experimentally, pentane has been observed to form a thin 
wetting film on a water substrate (15 ) . Films of octane on water are marginally 
stable. A theoretical study of hydrocarbon adsorption on water surfaces that 
describes and confirms this phenomenon has been done using Lifshitz's theory (16) . 
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The values of £(i£j) for ideal metals (smooth and clean surfaces free from 
oxides and contaminates) are also presented in Pig. (XXI. 1). Since the complex 
part of the dielectric constant, 2n k, is large for metals, e(i£) also has a 
relatively large value. In practical heat transfer equipment these values 
cannot be reached since non-idealities (oxides and adsorbed contaminates) are 
usually present. However they do represent an upper limit for E(i^)* Experi- 
mentally, the important value is the measured value of 2nk = £" which character- 
izes the surface being used. These results also emphasize the fact that metals 
are high energy surfaces which are easily wetted and easily contaminated. 

The dielectric properties of polytetraf luoroethylene (Teflon) fall above 

but close to those of pentane. Pentane wets Teflon since its surface tension 
ci 2 

(y = 0.016 J/m ) is less than the critical surface tension for the wetting 

2 

of Teflon (0.0185) J/m ). From a more fundamental point of view, pentane wets 
Teflon because E(i£j) for pentane is less than e(i£) for Teflon. High energy 
diamond and graphite surfaces have also been studied extensively and the resulting 
data for E(i£) are presented in Fig. (III.l). 

Although glass is usually considered to be a high surface energy material, 
the van der Waals interaction constant for g. ass turns out to be close to those 
of typical polymer surfaces (17) . Therefore, the dielectric properties are 
not presented in Fig. (XII. 1). This is confirmed experimentally since the 
wettability of untreated glass surfaces is poor (18). However, the surface 

f 

characteristics and wettability of glass surfaces can be changed by various 
treatments (18). As a result, the surface history (cleaning and preparation) 
of the glass surface has a considerable influence on the experimentally deter- 
mined contact angle and many apparent inconsistencies are present in the litera- 


ture . 



materials. 
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III. 5 Theoretical Values of the Hamaker Constant, A 

— s£v 

Churaev used Eq. (III. 4) to calculate the Hamaker constant for various 

systems (22). These results are presented in Table (III.l). Since the dielectric 

constants for metals are larger than that for water which, in turn, is larger 

than that for pentane , the absolute value of the Hamaker constants for the 

pentane - metal systems are larger than that for the pentane - water system. 

~ -21 

For comparison the Hamaker constant for the He - CaF^ system is A s £ v ~“3.8 x 10 

J. (20). The lower values of the dielectric constants for Helium and CaF^ 

result in a relatively low value of A . . The refractive index of CaF„ is 

s£v 2 

approximate ly n ~ 1.4. 

The approximate value of the Hamaker constant, A ^ v , can be obtained 

using 


A s£v “ A U “ A s£ 


(III. 19) 


0.5 


A good approximation for the Hamaker constant, A^, is (A A^) * (23). 


Therefore we can calculate from the more readily available values of 

A ss and Ap ff . Calculated values for obtained using Eq. (III. 19) are 


"££ 


s£v 


presented in Table (III. 2). The values for A^ were obtained using 


A ££ = 24Tr d Y 


(III. 20) 


Israelachvili (24) successfully compared the surface tension of a number of 
saturated hydrocarbons calculated on the Basis of Eq. (Ill .20) with experi- 
mental data. The value for gold was obtained from Ref . (23) and that for Teflon 

from Ref. (17). in Table (III. 3) a list of values for A for other solids from 

ss 

these references is given. Graphs of Eq. (III. 19)' are presented in Fig. (III. 2) 

-19 -19 

for various solids: Au, A = 4.5 x 10 J ; A£_ 0^ , A = 1.5 x 10 J; 

ss 2 3 ss 

-19 

polystyrene, A ss = 0.65 x 10 J. The values for A gs were obtained in Ref. (23). 


Due to limitations in the availability and interpretation of the optical data 


there is a degree of uncertainty in the values of these constants. For example, 

-19 

the curve for A = 1.05 x 10 J represents the value for Teflon obtained 

by Vassi"' ieff and Ivanov (17). Their value for Polystyrene is A £s = 1.11 x 10 ^ 

-19 

which is considerably different from the value A = 0.652 x 10 ' J referenced 

ss 

above (23). As a result of these inconsistencies, the theoretical results based 
on these initial values of the Hamaker constant are obviously approximate. 
However, the relative values calculated using the same set of assumptions 
are very instructive. The refinement of the procedures and the availability 
of more extensive optical data will lead to more consistent results in the 
future. The potential use in heat transfer will be an added stimulus for 
obtaining more extensive optical data. 


Table (III.l) 


Calculated Values of A ^ for Wetting Films (22) 



Liquid, substrate 

-w- 20 

Liquid, substrate 

20 

- A s£v * 10 



J 


J 


Water on quartz 

1.12 

Pentane on water 

0.74 

t 

Octane on quartz 

1.67 

Hexane on water 

0.63 

It 

Decane on quartz 

1.57 

Decane on water 

0.34 

:{ 

Tetradecane on quarts 

1.34 

Octane on steel 

< 

19.2 

}; 

Benzene on quartz 

0.91 

Decane on steel 

18.9 

r? ■ 

C C&a on quartz 

0.56 

Decane on gold 

23.1 


Water on gold 


11 • 2 



Table (III .2) 

Values of A „ Calculated Using Eq. (III. 19) 


1rt 20 
A *10 
ss 

J 

v 10 ' 0 

J 

,„20 

- A stv- 10 

J 

Gold : 45 

Methane : 3.95 

9.4 

Gold : 45 

Pentane : 5.04 

10.0 

Gold : 45 

Hexane : 5.48 

10.2 

Gold : 45 

Heptane : 5.82 

10.4 

Gold : 45 

Octane : 6.08 

10.5 

Gold ; 45 

n*Dodecane : 6.75 

10.7 

Gold : 45 

Nitrogen : 2.49 

8.1 

Teflon 10.5 

Methane : 3.95 

2.5 

Teflon : 10.5 

Octane : 6.08 

1.9 

Teflon : 10.5 

Nitrogen : 2.49 

2.6 


Table (III. 3) 

Calculated Values of A. . from Ref. (17) 


Material 

A. . *10 2 ° -J 
n 

Material 

A. . *10 20 J 
ii 

Decane 

6.7 

Poly propylene 

9.05 

Diamond 

41.4 

Teflon 

10.5 

Hydrocarbon II 

7.7 

Quartz 

8.65 

Polyethylene 

Polystyrene 

10.5 
11 1 

Water 

6.05 
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FIGURE (III. 2) Calculated value of the Hamaker constant for a thin film of 
liquid on a solid substrate (Eq. III. 19). 




1 
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III. 6 Theoretical Values of the Interline Heat Flow Nuraher, A h „ /V 

fg 

The product of the surface transfer coefficient, A/V , and the heat of 

vaporization, h , gives a measure of the interline heat sink capability of 

the liquid - solid system, which we call the Interline Heat Flow Number, h_ A/V. 

fg 

For the constant heat flux model, the average heat flux over the meniscus thick- 
ness range r) > 1 is obtained using (Eq. II* 15c) 

A h. 


q = 


fg 


x 


V 


2 a r) 


(III. 20) 


The lower limit for the length of the thin film, x, of thickness range, r) > 1, 
depends on meniscus stability and is beyond the scope of this study. Assuming 
bulk values for the heat of vaporization and the kinematic viscosity, the inter- 
line heat flow number can be calculated using the values for A = -A s ^ v /6tt 
previously obtained. Theoretical values for the interline heat flow number for 
a few representative systems are presented in Table (III. 4). The gold-alkane 
system and the Teflon -alkane system are two ideal systems for comparison since 
the solids have been extensively studied and the a^Kanes are simple non-polar 
fluids. More importantly, these systems represent a broad range since gold 
is a very high surface energy material and Teflon is a very low surface energy 

material. The value of the interline heat flow number for methane-gold is 
-9 -9 

9.76 x 10 W versus 2.6 x 10 W for methane-Teflon . Increasing the molecular 

f 

— -9 

weight of the fluid gives h^ A/V = 2.62 x 10 W for octane-gold versus 

-9 

0.26 .c 10 W for octane-Teflon. This decrease occurs because the ratio 

h_ /V decreases faster than A increases. For all practical purposes these’ 
fg 

numbers bracket the range for simple fluids. The inconsistencies present in 
the table values for the interline heat flow number result from the inconsis- 
tencies in the values of A in the literature in the different references. For 
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Table (III. 4) 


Theoretical Values of Interline Heat Flow Number, h„ A V 

fg 


Liquid - Substrate 

(h 

fg 

- -1 9 
A V ) *10 

Watts 

A 

Source 

Methane - Gold 

9.76 

(111.4K) (1) 

Table 

(III. 2) 

Pentane - Gold 

5.26 

(293K) 

Table 

(III. 2) 

Hexane - Gold 

4.22 

(293K) 

Table 

(III. 2) 

Heptane - Gold 

3.33 

(293K) 

Table 

(III. 2) 

Octane - Gold 

2.62 

(293K) 

Table 

(III. 2) 

Do decane - Gold 

1.12 

(293K) 

Table 

(III. 2) 

Nitrogen - Gold 

4.17 

(77K) 

Table 

(III. 2) 

Decane - Steel 

2.90 

(293K) 

Ref. 

(22) 

Decane - Steel 

4.14 

(447K) (2) 

Ref. 

(22) 

Decane - Gold 

3.50 

(293K) 

Ref. 

(22) 

Methane - Teflon 

2.60 

(111.4K) 

Table 

(III. 2) 

Octane - Teflon 

0.26 

(293K) 

Table 

(III. 2) 

Nitrogen - Teflon 

1.53 

(77K) 

Table 

(III. 2) 

C Ci. - Quartz 

0.105 

(293K) 

Ref. 

(22) 

4 

C Cl, - Quartz 
4 

0.497 

(2) 

(350KV 

Ref. 

(22) 

Octane - Quartz 

0.418 

(293K) 

Ref. 

(22) 


(1) and V are evaluated at this temperature. 

fg 

Ref. (22) values for A were adjusted for temperature using . 
Eqs . (III. 19 and III. 20). 


( 2 ) 
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IV. The Theoretical Interline Heat Sink, Q 

The above results can be used, to estimate various characteristics of the 
interline evaporation process. The theoretical interline heat sink for a unit 
length of interline, Q, is obtained using 

Q = q x = [2 In n]°' 5 [q] 0 ' 5 [h fg A/V]° -5 (IV. 1) 

The length of the interline region of thickness range q > 1, x, is 

x = [2 in n]°' 5 [q]“°‘ 5 [h fg A/v] 0,5 (11.13) 

These equations are presented in Figs. (IV. 1 and IV. 2) for n = 10. In Fig. (IV. 1) 
some of the values of Ah. /V listed in Table (III. 4) are marked on the theoret- 

fg 

ical curve. Obviously an increase in the interline heat flow number gives an 

increase in the theoretical heat sink, Q, and the length cf the interline region, 

x. Two of the systems, decane-steel and carbon tetrachloride-glass are given 

in more detail in Fig. (IV. 2). Curve (la) represents the decane-steel system 

at 293K and Curve (lb) represents the system at 447K. Curve (2) represents the 

carbon tetrachloride system at 293K. The results predict that the interline 

heat sink, Q, for the decane-steel system is approximately three times larger 

than that for the carbon tetrachloride system at the same heat flux and 
5 -2 

P = 1.01 x 10 N*m . This results from the increased length of the evaporating 
v 

thin film, x, of thickness 1 < n 10, which is a function of the interline heat 
flow number. The increase of Q with temperature results from the substantial 
decrease in the kinematic viscosity in this temperature range. This increase 
will not continue since the dispersion constant will decrease substantially at 
higher temperatures. Fcr the carbon tetrachloride-quartz system, both the 
dispersion constant and the reciprocal of the kinematic viscosity increase with 
temperature at this temperature level (note the effect of the relative size of 


i 


t 



Ul 

KJ 


FIGURE (IV. 1) Graphical representation of Eq. (IV. 1) and Eq. (IV. 2). 
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FIGIJKE (IV. 2) Heat sink capability, Q, and interline length, x, for n - 10; (la) decane-steel at 293K» 
(Jh) decane-steul at 44/K; (2a) CCV - quartz at 293K; (2h) CCt.. - quartz at J50K. 
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A on A ^ in Fig. (III. 2}) 


For a more extensive numerical example we compare the heat transfer 


characteristics of an evaporating thin film of carbon tetrachloride on quartz 


with those of decane on steel in Table (III. 5). The size of the ideal heat 


flux is approximately fixed by the vapor pressure and temperature difference 


for all practical purposes. Therefore, it is desirable to compare the two 


systems at the same vapor pressure level for which we use the atmospheric 


5 -2 

pressure of 1.01 x 10 N*m . Although this is desirable from a heat transfer 


point of view, it does add an additional uncertainty since the extrapolated 


values of A h V presented in Table (III. 4) are used. The interline tempera- 


-2 

ture difference of 10 K is used to fix the heat flux level at approximately 


q = 16,000 W*m . Due to the larger value of A, the interline thickness, 


■6 q , ar i the interline length, x, of the decane system are larger than those of 


the carbon tetrachloride system. This gives the decane system a higher heat sink 


capability, Q. The larger values of the interline heat flow number and the 


thermal conductivity of the solid combine to give the decane-steel system a 


much larger effective overall heat transfer coefficient, For comparison. 


the heat transfer characteristics are also given for T = 293K. Although it is 


not possible to predict the maximum stable heat flux at the present time, the 


larger values of A and for the decane-steel system should also give a larger 


maximum stable heat flux. 


Since the constant heat flux model, Eg. (II. 2) , does not include a unique 


value of the liquid-solid interfacial temperature difference at the interline, 


(T„ - T )„, the same numerical results can be obtained using Eas. (11.13) and 

£v v 0 


(11.15) with a different set of interline conditions. For example, keeping r) 


and a unchanged, the same numerical results obtained in Table (III. 5) for x and Q 


with CC£.-cmartz can be obtained using Y = 0.05 and (T 0 - T ) = 0.1. In this 

4 * £v v 0 


example, the average interfacial temperature difference which is fixed by the 
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Table (III. 5 ) 

Comparison of Carbon Tetrachloride-Quartz System 
with the Decane-Steel System 



» 10 
o *10 , m 

o 


„ -10 -2 
G* 10 , m 


, id -4 -2 -1 

h n *10 , M*m *K 

£v 


— 22 
A* 10 , N*m 


h. A 

-2- *lo 10 , W 

V 


-4 -2 

P • 10 , N*m 

v 


—id , -4 „ ~2 

q *10 , W*m 


— -4 -2 

q* 10 , W*m 


3 -1 

Q* 10 , W-m 


x*10 , m 


4 -] —l 

U *10 , Wm *K 
eff 


{T - T ) , *K 
s v 


(T *v ' V ' K 

o 


C C£^ - Quartz 
T - 293K 

V 

C 10 H 22 " SteSl 
T = 293K 

C 10 H 22 " Steel 
T v = 447. 2K 

29 

103 

121 

7,000 

6.4 

800 

74 

1.85 

330 

3 

100 

66 

1.05 

29 

41.4 

1.19 

0.012 

10.1 

0.74 

0.0186 

3.3 

0.37 

0.0093 

1.65 

0.5 

0.5 

0.5 

10 

10 

10 

1. 34 

1.12 

17.5 

J .36 

12 

1.05 

3.10 

2,040 

179.8 

io ~ 3 

io ~ 3 

io ~ 3 

=2.85 

=0.0036 

=0.65 

0.01 

0.01 

0.01 



6,560 


212.7 
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heat flux, (T^ - T ) = g /h^, is also unchanged. However, the film thickness 

0 - o — 

range becomes 14 <_ o _< 140A instead of 29 <_ o <_ 290 i. Therefore , y is a measure 
of the tempera tue change along the solid surface in the immediate vicinity of 
the interline which is not included in the constant heat flux model. The con- 
stant temperature model demonstrates that a large percent of this temperature 
change occurs in the thickness range 1 <_ p <_ 2. Which, in turn, indicates that 
the average heat flux, q, fixes the average value of (T^ - T^) at approximately 
its ideal value except in the region 1 <_ n <_ 2. in this region the local tempera- 
ture gradient with (T„ - T ) > (T„ - T ) fixes 6 . 

3 £v v 0 x,v v 0 


vj 


I 1 ! I 
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V. Discussion 

A theoretical procedure to determine the interline heat sink capability of 
an evaporating meniscus using the macroscopic optical properties and the thermo- 
physical properties of the system has been outlined above. This procedure signi- 
ficantly increases our understanding of the interline and shows how the inter- 
line transport processes are controlled by the dispersion forces resulting from 
the interaction of the solid and liquid. The procedure can be used to predict 
the relative size of the interline heat sink and, thereby, enhance the develop- 
ment of heat transfer equipment which use tile evaporating meniscus and/or the 
interline region. 

Due to the basic nature of the equations, their current level of development 
and the limited availability of optical data, the use of ideal systems 
(e.g. neglecting, surface roughness) has been emphasized. However, this does not 
detract from the significance of the results since ideal systems give an order 
of magnitude approximation and enhanced understanding of real systems , and 
reveal the basic mechanisms involved in the transport processes . 

The constant heat flux model explicitly demonstrated the importance of 
the interline heat flow number. Using theoretical values of the interline heat 
flow number, the relative size of the interline heat sink and the interline 
heat transfer coefficient can be determined. The absolute upper limit of the 

r 

interline heat sink cannot be determined using only these procedures because 

the important question of meniscus stability, which must be evaluated experi- 
mentally, has not been addressed herein. However, these concepts do describe 
why the interline heat sink of a system like decane-steel should be greater 
than that of a system like carbon tetrachloride-quartz. A procedure is now 
available whereby various systems can be theoretically compared. 

By necessity, there are aspects of the above theoretical results which 
need experimental evaluation, refinement and extension because of the 
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assumptions in the models. These results indicate that this can be accomplished 
by measuring the meniscus stability and the interline length, x, as a function 
of heat flux and the interline heat flow number. 
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Appendix A 

Determination of Values of £(ig) Used in Fig. (III.1) 

There are various procedures to obtain spectral graphs of £ (i£) from 
the limited experimental data. One of the most useful makes use of the 
following empirical formula (19) 

fjiffrr - * i- b5 > w- 1 ’ 

This equation was used to obtain the curves for graphite , diamond and Teflon 

presented in Fig. (III.l). The constants and references are given in Table (A.l) . 

The curve for pentane was obtained from Ref. (16) . The dielectric constant 

2 

is essentially constant (£ = n ) from zero frequency through to an ultraviolet 
frequency, to , which they choose to be the first ionization potential. For 
frequencies less than U) they used the following equation: 

2 

e(i£) = 1 + - - — V (A- 2) 

■ 1 + (C/UJ ) 

uv 

For frequencies higher than the plasma frequency, w t they used Eq. (A. 3) 

JP 

with straight line interpolation in between 

e(i5) = l + w 2 /^ 2 (a-3) 

p 

r 

The curve for Helium was obtained from Ref. (20) in which a similar procedure 
was used. 

The data for nitrogen, oxygen, methane and carbon tetraf luoride was 
obtained using 


(A-4) 


r 


and standard values for the index of refraction. The values for gold and 
silver were obtained from Ref. (21) . 


Table (fl.l) 

Constants and References for Eq. (A.l) 


Substance 

Graphite 

Diamond 

Polyethylene 

Polypropylene 

Teflon 

Polystyrene 


a 

b'lO'^, s/rad 

Rei 

0.948 

3.9 

19 

0.705 

1.5 

17 

0.420 

2.1 

17 

0.390 

2.1 

17 

0.317 

1.2 

17 


0.433 


17 
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